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ABSTRACT 

A self-similar formalism for the study of the gravitational collapse of molecular 
gas provides an important theoretical framework from which to explore the dy- 
namics of star formation. Motivated by the presence of elongated and filamentary 
structures observed in giant molecular clouds, we build upon the existing body 
of work on cylindrical self-similar collapse flows by including dynamic equations 
of state that are different from the effective equation of state that produces the 
initial density distribution. We focus primarily on the collapse of initial states 
for which the gas is at rest and everywhere overdense from its corresponding 
hydrostatic equilibrium profile by a factor A, and apply our results toward the 
analysis of star formation within dense, elongated molecular cores. An important 
aspect of this work is the determination of the mass infall rates over a range of 
the parameters which define the overall state of the gas - the overdensity pa- 
rameter A, the index F of the static equation of state, and the index 7 of the 
dynamic equation of state. While most of the parameter space explored in this 
work leads to solutions for which the underlying equations do not become sin- 
gular, we do include a discussion on how to treat cases for which solutions pass 
smoothly through the singular surface. In addition, we also present a different 
class of collapse solutions for the special case 7 = 1. 

Subject headings: hydrodynamics - stars: formation 
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1. INTRODUCTION 

Star formation in our galaxy occurs primarily within giant molecular clouds (GMC) - 
highly nonuniform complexes of molecular gas containing a total mass of ~ 10^ Mr^ within a 
radius of ~ 20 pc. These complexes have hierarchical structure that can be characterized in 
terms of clumps and dense cores surrounded by an interclump gas of density ~ 5-25 cm~^. 
Clumps have characteristic densities of ~ 10^ cm~^ and radii ranging between 0.2-2 pc, 
the largest of which are comprised of as many as ~ 1000 small {R ~ 0.1 — 0.2 pc), dense 
(~ 10''-10^ cm""') cores whose mass function has been measured to range from ~ 1 — IOOMq 
(with a peak ~ IOMq) by Jijina ct al. (1999), and more recently, from ~ 0.2 — 2OM0 (with a 
characteristic mass of ~ 2M0) by Lada et al. (2008). It is the gravitational collapse of these 
cores within a clump that results in the formation of stars. The dynamics of core collapse is 
therefore a fundamental component of the star formation process. 

The study of self-similar collapse flows has provided an important cornerstone of the 
current theory of star formation (e.g., Shu et al. 1987), with numerous works appearing 
in the literature. The original self-similar collapse calculations (Larson 1969ab; Penston 
1969ab; Shu 1977) considered isothermal spherical flows. Since then, many generalizations 
of the collapse have been made. The leading order effects of rotation have been studied, 
both for the inner pressure free region (Ulrich 1976; Cassen & Moosman 1981) and for the 
entire core (Terebey, Shu, & Cassen 1984). The leading order effects of magnetic fields have 
also been included (GaUi & Shu 1993ab; Li & Shu 1996, 1997). More recently, the collapse 
of magnetized singular isothermal toroids has been studied (Allen, Shu, & Li 2003; Allen, 
Li, & Shu 2003). 

While much of the focus has been given to spherically symmetric flows, fllamentary 
and elongated structures in molecular clouds arc commonly observed (e.g., Houlahan & 
Scalo 1992; Harjunpaa et al. 1999; Jijina 1999), and based on the results of numerical 
simulations (e.g., Curry 2000; Jappsen et al 2005), appear to be an important aspect of the 
star formation process. A complete understanding of molecular cloud dynamics thus requires 
the inclusion of cylindrical geometries. Indeed, several authors have applied self-similar 
techniques toward the study of how cylindrical structures collapse (Inutsuka & Miyama 
1992; Kawachi & Hanawa 1998; Hennebelle 2003; Tilley & Pudritz 2003; Shadmehri 2005). 
A central element of these studies involves the form of the equation of state used. While 
an isothermal equation of state clearly provides a reasonable and important starting point 
for the general study of core collapse, observational evidence points to softer equations of 
state. Specifically, non-thermal linewidths Aw in molecular cloud cores show a correlation 
with density of the form (Af)^ ~ p~^, with (3 > (e.g., Larson 1981; Jijina et al. 1999). 
If one interprets the linewidth Av as the effective transport speed in the medium, then the 
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corresponding effective equation of state has the form P , where T — 1 — (3 and hence 

< r < 1. 

The self-similar collapse of filamentary structures with polytropic equations of state has 
been explored by Kawachi & Hanawa (1998) and Shadmehri (2005), with the latter work also 
including the effects of magnetic fields. These works, however, considered static equations 
of state which imposed global constraints on the pressure and density. It is likely that the 
physical processes which govern a gas will change as that gas undergoes gravitational collapse. 
We therefore extend the analyses of previous works on cylindrical flows by considering a 
dynamic equation of state (Fatuzzo, Adams, & Myers 2004). Specifically, wc consider the 
case in which the dynamic equation of state for the collapsing gas is different than the 
effective (static) equation of state that produces the initial equilibrium configuration. Here, 
the static equation of state (as set by F) refers to the pressure law that enforces the initial 
(pre-coUapse) configuration for the gas, whereas the dynamic equation of state (as set by 7) 
refers to the pressure law that describes how the thermodynamic variables of the gas change 
as the material is compressed during collapse. This process is governed by the entropy 
evolution equation (introduced in §2.1), so that specific entropy is conserved along a given 
streamline. However, since the physics that determines the density profiles of the pre-collapse 
states can be different from the physics that governs the thermodynamics of the collapse ffow, 
we allow 7 7^ F. Our analysis therefore allows for a more robust model from which to gain 
insight into the dynamics of filamentary/elongated molecular structures. 

An important aspect of this work is its focus on the mass infall rate M . Collapse fiows 
are often self-similar and have no characteristic mass scale. However, the mass infall rate is 
one of the most important physical quantities in the star formation problem, and determines, 
in part, the total system luminosity and the total column density of the inf ailing envelope. 
These quantities, in turn, largely account for the spectral appearance of protostellar objects 
(e.g., Adams, Lada, & Shu 1987; Adams 1990) since most of the luminosity is derived 
from material falling through the gravitational potential well of the star. Although the 
circumstellar disk stores some of the energy in rotational motion, the system luminosity is 
(usually) a substantial fraction of the total available luminosity 

-rt* 

where Mg is the total mass of the star/disk system, M is the mass infall rate, and i?* is the 
stellar radius. The stellar radius, which helps determine the depth of the potential well, is 
itself a function of the mass infall rate (Stabler, Shu, & Taam 1980). 

The paper is organized as follows. We formulate the collapse problem via self-similar 
methods for the general case of a collapsing cyhndrical cloud of gas in §2. In §3, we determine 
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the range of parameter space that yields collapse solutions by considering the limiting cases 
v/x ^ and x/v ^ 0, and illustrate how the ensuing collapse is affected by the initial state 
of the gas. Guided by these results, we explore in §4 the collapse of initial states that are out 
of exact hydrostatic equilibrium by being overdense, and apply our results to the collapse 
of elongated cores in §5. We consider cases for which solutions go smoothly through the 
singular surface in §6, including the 7 = 1 case which leads to a different class of collapse 
solutions. We present our conclusions in §7. 



2. FORMULATION OF THE COLLAPSE PROBLEM 

2.1. Basic Governing Equations 

We adopt a cylindrical coordinate system described by the variables r, z and 0. To keep 
the problem tractable, we assume dependence only in r (i.e., the cylindrical structure is infi- 
nite in length) , with the self-gravitating gas described locally by its density p(r, t) , pressure 
P{r,t) and (radial) velocity u{r,t), and globally by the mass per unit length ML{r,t) con- 
tained within a radius r. The gravitational collapse of this fluid is governed by conservation 
of mass 

-^ + u^ = and ^ = 2.rp, (2) 
or equivalently, by the equation of continuity. 



dp 1 d 

and the force equation 



du du IdP 2GMl 

^ + '"^ = — ^ ■ (^) 

at or p or r 

To complete the description required for the evolution of the gas to be solved, the 
pressure must be specified through a choice of the equation of state. For example, the 
relation P = s^p is used to describe an isothermal gas (where s is the sound speed), whereas 
a polytropic equation of state P = Kp"' allows for a more general treatment of the problem. 
Of course, adopting an equation of state represents a simplification to the real system as 
it embodies the numerous physical processes which govern the true state of the gas - and 
how they add/remove energy to/from the gas - into one simple relation. It is quite likely 
that these processes will change significantly during the collapse, so that the equation of 
state which governs the gas during the initial stages of collapse will almost certainly have a 
different form from that which governs the collapse during the later stages of the collapse. 
We therefore introduce a dynamic polytropic equation of state in our formalism which allows 
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the relation between pressure and density to evolve during the collapse (e.g., Fatuzzo, Adams 
& Myers 2004). Specifically, we assume that entropy is conserved along a streamline, and 
use the conservation of entropy equation 

to follow the evolution of the pressure. We thus refer to 7 as the index of the dynamic 
equation of state. It is important to note that equation describes how a given parcel 
of gas changes its thermodynamic variables along a streamline, therefore allowing for an 
equation of state which evolves during the collapse. In contrast, relating the pressure to the 
density through a fixed (static) equation of state (e.g., P = ICp^) implies a global constraint 
on those variables. 



2.2. The Similarity Transformation 

As shown in §2.1, the cylindrical collapse problem is represented mathematically by a 
set of coupled partial differential equations in time t and radial position r. In this section, we 
find a similarity transformation that reduces this set of PDE's to a set of ordinary differential 
equations in a new similarity variable x which we define below. In particular, we look for a 
similarity transformation of the general form 

X = Afr , p = Bt''a{x) , Ml = Cfm{x) , 

u = Dt%{x) , and P = Efp{x) . (6) 

Here, both the coefficients (A, B, C, D, E) and the indices (a, 6, c, d, e) are constants. The 
reduced fluid fields (a, m, v, p) are dimensionless functions of the (single) dimensionless 
similarity variable x. In this work, the time benchmark t = corresponds to the instant of 
the onset of collapse. 

The general similarity transformation calculation leads to four equations to specify the 
five indices a, b, c, d, e. We leave the constant a arbitrary for the moment and write the rest 
of the variables in terms of its value, i.e., 

a = a, b = —2 , c = — (2a + 2), 

d = -{a + l), and e = -2(a + 2) . (7) 
Similarly, for the coefficients we obtain 



A = A, B = {2ttG)-\ C = {A^G)-\ 
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D = A-\ and E = {2nGA^)-^ , (8) 

where G is the gravitational constant. We thus obtain reduced equations of motion in the 
form 

{ax + v)— = {2a + 2)m, (9) 
dm 

— — = xa , (10) 

dx 

, ,1 da dv / v\ ,^ ^ , 

(^ax + v)-— + — =(2--) , (11) 
a dx dx V X/ 

^dv 1 dp 2m , ^ , ^ 

{ax + v)— + — f = + {a + l)v, 12 

dx adx X 

iax + v)^\og[p/a^] = 2{2 + a--f). (13) 

Note that our similarity transformation is not unique - one can always rescale the coefficients 
{A, B, C, D, E} by a set of dimensionless numbers and obtain new equations of motion with 
different numerical coefficients. 

The first two equations of motion can be immediately combined to obtain an expression 
for the reduced mass m{x), i.e., 

(ax + v) 

m = - -xa . 14 

(2a + 2) ^ ^ 

Likewise, multiplying equation IQ by the constant 

^ = ^r^^2 + a-7). (15) 
2a + 2 



and subtracting the product from equation (]T3l) yields a differential equation which can be 
integrated to obtain an expression for the reduced pressure 



p= Cia^m" = Cia<^'i 



where Ci is a positive integration constant. 



{ax + v) 
(2a + 2) ^_ 



1 



(16) 



Given the solutions for the reduced pressure p{x) and reduced mass m(x), Eqs. (fTTl) 
and (IT^ are the relevant equations of motion to determine the remaining unknown functions 
a{x) and v{x). Using Cramer's rule, we derive an equivalent set of equations 

da A ^ dv V . 

—r- = and -;- = ^ , (17) 

dx V dx V ^ ' 
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where 

^^{ax + vl_jp (18) 

V^-^^^- '^^'+r^^ +K + ^)a + »)--^(2-^), (19) 
1 + a a \ X/ 

and 

/ ^ /'^ v\ a(ax + v) 2p (2 + a — 7) , 

^ ^ \ xJ 1 + a a{ax + v) ^ ' ^ ' 

The physical state of the gas is defined through a choice of the parameters (a, 7, Ci), and 
its collapse is described by the solutions to Eqs. (14) - (20) for the specified set of reduced 
field variables Vi — v{xi) and ctj — a{xi). 

It is mathematically possible to consider "complete" solutions that span the entire 
available range —00 < x < 00 (first obtained by Hunter 1977; see also Whitworth & Summers 
1985), and as such, span both negative and positive times. However, it is rather unlikely 
that molecular clouds will evolve toward their centrally condensed initial configurations in a 
self-similar manner subject to (only) the physics included in these equations of motion. For 
example, before the onset of collapse (for t < 0), molecular clouds may evolve through the 
processes of ambipolar diffusion (at least for small mass scales), shocks, turbulent dissipation, 
cooling fiows, condensation instabilities, and cloud-cloud collisions. In addition, the cloud 
will most likely initiate collapse before a completely self-similar equilibrium state has been 
attained; the collapse will only become self-similar asymptotically in time (i.e., the self- 
similar collapse solutions of this paper are intermediate asymptotic solutions to the realistic 
problem of the collapse of a finite cloud with finite central density). As such, the self-similar 
solutions of this paper for the protostellar collapse phase (t > 0) cannot (in general) be 
extended to the prc-stellar phase (t < 0), as they would likely encounter a critical point and 
become singular. We therefore limit this discussion to solutions with < x < 00, sometimes 
called "semi-complete" solutions. 



3. GENERAL SOLUTIONS TO THE COLLAPSE PROBLEM 

In order to determine what set of parameters (a, 7, Ci) yield solutions that describe the 
collapse of a molecular cloud filament, we obtain analytic solutions to the equations of motion 
in the limit that t — > 0. As these solutions describe the early stages of the collapse, the gas 
velocity u must be negative and bounded. For a positive value of a, this constraint requires 
that u — Dt~'^"''^^^v{x) not be singular at t — 0, which in turn requires that v{x) oc oc t^", 
where P > {a + l)/a. Clearly then, the ratio — > as i — > for solutions relevant to our 
discussion. 
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In this limit, Eq. (9) reduces to a form that can be easily integrated, and along with 
Eqs. (14) and (16), leads to analytic solutions for the reduced density, mass, and pressure 
of the form 

a = Ax2/", (21) 



2a + 2 



and 



p = Ci 



1 1 



a 



2+a 



(22) 



(23) 



_2a + 2_ 

where A is a positive constant. The reduced velocity is governed by the limiting form of Eq. 
(12) 

dv 

-(a + i; 

X 

where 



a— -(a + l)- = yoa:^^ , 



K = -2A 



1 + Ci 



2 + a\ /2a + 2 



_2a + 2 

Eq. (24) yields a power-law form solution 



A" 



(7-l)/(l+a)' 



(24) 
(25) 



(26) 



indicating that collapse solutions can be found for a > (since Vq < 0). However, the 
real density becomes time-independant and scales with radius as p oc r^/" in the limit that 
x/v — > 0. Since the density approaches zero at small radii, solutions for a > cannot 
represent the collapse of filamentary like structures (whose initial density profiles are peaked 
at r = 0), thereby ruling out this range of parameter space in our work. 

Interestingly, the ratio v/x also approaches zero as i — > for the case that a < (since 
X oo). The reduced density, mass and pressure given by Eqs. (21) - (23) therefore also 
describe the state of the gas at the onset of collapse for this case. However, the real density 
(which still scales as p oc r^/") now becomes singular as r — 0. Of course, the equations 
of motion presented in this work are simply mathematical idealizations to the real, physical 
problem, for which filaments have a finite density core at small radii r < tq and a large but 
finite outer boundary Tout- This "real" system is expected to follow the self-similar solutions 
for intermediate length-scales. Indeed, previous numerical work (Foster & Chevaher 1993) 
indicates that the collapse of an isothermal core approaches the expected self-similar form 
when the core has routine > 20. Likewise, Fatuzzo, Adams & Myers (2004) found similar 
results, with cores that have initial inward velocities more readily approaching the self-similar 
collapse forms. 
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For the case that a is negative, additional constraints are provided by the requirement 
that the reduced mass m remain positive, and that the reduced pressure p increases as a 
increases. Together, these physical conditions require —2 < a < — 1, as can be clearly seen 
by the forms of Eqs. (22) and (23). In addition, collapse solutions (for which Vq < 0) require 
that 



A > A 



crit 



1/a 



(2 + a) 



(a+l)/(a-a7) 



(27) 



_2a + 2_ 

We note that Xcrit is not defined if 7 = 1. Indeed, the 7 = 1 case yields a unique class of 
collapse solutions, and will be considered separately in §6.2. 

Since the real pressure is defined in terms of the constant A through the similarity 
transforms, we can set p — 0;"+^ without loss of generality. That is, for a real pressure 



setting 



A = [X:(27rG)i-^ 



i-ri-1/2 



and r = (2 + a) yields the desired form for the reduced equation of state {p 
sets the integration constant 

„ "2a + 2"" 
(^1 



A"^ 



and, in turn. 



(2a + 2)(2 + a) 



l/a 



(28) 

(29) 

0;"+^), and 

(30) 
(31) 



To obtain full solutions, one must numerically integrate Eqs. (17) - (20) from a specified 
initial set of values ctj = a{xi) and Vi — v{xi). We note that one cannot set initial conditions 
which correspond to the onset of collapse since x — > 00 as t — > 0. As such, initial conditions 
are set through an arbitrary choice of >> 1 along with a corresponding value of ctj for 
which A > Xcrit s-iid a small but negative value of f j. Representative solutions arc illustrated 
in Figs. (1) and (2) for the case that a = —1.25 (F = 0.75) and 7 = 0.25 (with Ci given 
by Eq. [30]), obtained by integrating inward from rcj = 5 x 10^ and a-j = 5.78 x 10"'^ (for 
which A = 1.5Acrjt)- The solid curves represent the solutions obtained by setting the value 
of Vi equal to its power-law counterpart, as defined by Eq. (26). The dashed lines represent 
solutions obtained by setting the value of Vi equal to 0.1, 0.3, 3 and 10 times the power-law 
value. It is clear from our results that the reduced density is not sensitive to the initial 
velocity, as all five cases yield values that differed by less than 0.1 percent. In addition, it is 
clear that the velocity evolves toward its power-law form as given by Eq. (26). 
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For completeness, we note that in the hmit that x/v ^ 0, the reduced mass approaches 
a constant value - m — > mo- In turn, the reduced pressure takes the form p oc a^. Taken 
together, the two limiting solutions presented here thus illustrate how the equation of state 
smoothly transforms from an initial index of F = a + 2 (at early times and/or large radii) 
to a dynamic index 7 (at late times and/or small radii). This result is clearly illustrated 
in Fig. 3, which plots p versus a for the solutions presented in Figs. (1) and (2). The five 
initial velocities yield results that are within 0.5 percent of each other. As can be easily seen 
from Fig. (3), the equation of state index evolves from a value of F = 0.75 (in the lower 
left corner of the figure) to a value 7 = 0.25 (in the upper right corner), with the transition 
occurring within a narrow region around x ~ 1. 

Finally, we note that the equations of motion which define our problem may well yield 
several different classes of solutions - both physical (e.g., wind solutions) and unphysical. 
However, solutions which pertain to the collapse of filamentary structures in the absence of 
shocks (but see §6) are all qualitatively similar, and therefore characterized by solutions such 
as those presented in Figs. (1) - (3). 



4. COLLAPSE OF OVERDENSE STATES 



The analysis presented in the previous section illustrates that collapse solutions of fil- 
amentary structures exist when — 2 < a < — 1 and A > Xcrit- It is easy to show that when 
A = Xcrit, Eqs. (21) - (23) yield self similar solutions which describe a gas in hydrostatic 
equilibrium (for which v — 0). These solutions are fully specified by the equation of state 
index F (which we use hereafter to specify the state of the gas instead of a), and reduce to 



aE 



mE = 



. (2 - 2r) 

(2-F)2 
f2-2r)r-i 



X2-T 



2-2r 
X 2-r 



(2 - rr 

and p = (assuming that A is set through the relation given by Eq. [29]). 



(32) 
(33) 



In turn, the real density and mass per unit length profiles are given by the time inde- 
pendent expressions 



Pe 





1 

2-r 




-2r)" 


1 

2-r _2 
7- 2-r 


2tiG_ 




[ (2 


_r)2_ 





and 



(2 - 2f; 



,r-i 



(2 - rr 



2-2r 
f 2-r 



(34) 



(35) 



In this section, we consider the collapse of gas structures which are initially at rest, but 
are out of hydrostatic equilibrium by being everywhere overdense by a factor A. As such, 
the initial gas (at time i = 0) is described by 



and 



0, 



(36) 



where a^; is the corresponding density profile if the gas were in hydrostatic equilibrium. 
Essentially, this class of solutions correspond to the solutions discussed in §3, where A = 
AAcrit, but with the velocity in the a; — > oo limit given by Eq. (26), as shown below. 



We note that hydrostatic equilibrium solutions of the form 



OLE 



and 



A 



a + 2 



■ X 



(2+2a)/a 



(37) 



exist when 7=1. However, since the reduced pressure takes the form p = Ciam, the 
corresponding balance between the pressure gradient and gravity, as described by Eq. (12), 

d 



Ci — (am) = -2 
ax 



am 



X 



(38) 



is then always maintained regardless of the value of A. The collapse scenarios considered in 
this section therefore cannot occur if the dynamic index 7 = 1, regardless of the value of 
the static index T. For completeness, we explore a different class of collapse solutions with 
7 = 1 in §6.2. 

With an initial (overdense) state defined in terms of F and A, as per Eqs. (32) and 
(36), and the dynamic equation of state defined in terms of 7, the ensuing collapse solutions 
can then be obtained by numerically integrating Eqs. (17). The adopted relation p = is 
maintained for the initial (overdense) state by setting 



(r-2)(7-r) 

A r^^r 



2-F 



7-r 
i-r 



(39) 



In practice, however, this initial state corresponds to a time tj = 0, and in turn, an 
initial value of Xi = 00. We therefore first obtain higher-order analytical expressions for the 
reduced density and velocity from those presented in §3 via a series expansion of the reduced 
equations in the limit Xq » 1 . Doing so yields the following values for the reduced density 
and velocity: 

AqXq- 



ao = Aasixo) 



-2 

. 2-r 



2-F 



(40) 



and 



-r 



2F3 
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where 



and 




/?o [1 - A-(^-^)] . 



2-r 



(42) 



(43) 



Full solutions can then be easily obtained by numerically integrating inward from xq. 

Wc present solutions for the reduced density and velocity in Figs. 4 - 7 for an overdensity 
parameter of A = 1.5 and several different sets of (7,r). The reduced density clearly 
exhibits the broken power law profile typical of self-similar collapse solutions. We note 
that the spectral indices which characterize the reduced density solutions are not sensitive 
to the choice of 7, and hence have similar forms when x « 1. In contrast, these solutions 
depend sensitively on the value of F, both in terms of shape when x » 1 and in overall 
normalization. Similar behavior is also clearly observed for the reduced velocity. 

As noted in the introduction, self-similar collapse flows have no characteristic mass 
scale. Instead, the collapse flow feeds material onto the central star/disk system at a well- 
deflned mass infall rate M. In these flows, the infalling material always approaches free-fall 
conditions on the inside (in the limit r — > 0) and the reduced mass determines the size of 
the infall rate. Specifically, m(x) — > mo as x ^ 0, and in turn, the mass accretion rate per 
unit length becomes 



in temporally increasing (decreasing) mass accretion rates. In contrast, the mass accretion 
rate for spherically symmetric flows is given by the expression 



and is constant in time for an initially isothermal gas (Fatuzzo, Adams & Myers 2004). 

We plot the value of mo as a function of F for A = 1.5 and three different values of 
7 in Fig. 8, and as a function of A for several sets of (F, 7) in Fig. 9 (the data point in 
this figure at A = 1.08 is discussed in §6.1). As expected, larger infall rates occur for initial 
states that are more overdense. In addition, stiffer static equations of state result in larger 
mass accretion rates. In contrast, the dynamic equation of state (as defined by 7) has little 
effect on the value of mo. 



Ml = 2(1 - F)/C {27rGy-^ Q-' mo t^'^^ . (44) 
Note that if F = 0.5, Ml is constant in time, whereas softer (stiffer) equations of state result 



M = (4 - 3F) /C=^/2 (A7rG)'^'-^y^ Q-' mo t'^'-^^ , 



(45) 



The insensitivity of the collapse dynamics to the value of the dynamic index 7 should 
not be surprising given the inside-out nature of the collapse (as shown explicitly in §5). 
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Specifically, the break at x fa 1 exhibited in the density and velocity profiles shown in Figs. 
4-7 denotes the boundary between gas that is in some part supported by pressure {x >> 1) 
and gas that has lost that support and is therefore approaching free-fall (x << 1). In terms 
of real quantities, this boundary occurs at a radius 

that moves outward in time. Since the speed at which this boundary moves is governed 
by the static equation of state, the collapse dynamics therefore depend sensitively on F. In 
contrast, the dynamic gas pressure has little effect on the collapse dynamics when x < 1 
(r < rs), and hence, on the overall collapse. 



5. APPLICATION TO ELONGATED CORES 

Dense star-forming cores in molecular clouds are on average shghtly elongated, with 
some cores exhibiting aspect ratios as large as 5 (Jijina et al. 1999). We apply our results 
to the evolution of fairly elongated cores, making the highly idealized assumption that their 
collapse can be reasonably approximated by our formalism. In reality, how such cores collapse 
likely falls within the limiting cases of cylindrical symmetry (explored here) and spherical 
symmetry (explored in Fatuzzo, Adams & Myers 2004). 

Interpreting the observed non-thermal linewidths Av in molecular cores as the effective 
transport speed in the medium, one finds 



Av 



dP 
'dp 



''^' = [/CF]^/'p-(l-^)/^ (47) 



for an assumed polytropic equation of state P = ICp^ . Adopting the fiducial values F = 0.5 
and Av = 1 km/s for a density of 5,000 cm^'^ (e.g., Larson 1981) then yields a value of 
K, — 2.6 g^/^ cm-*^/^ s~^, and in turn, density and velocity profiles given by 



p(r, t) = 2.4 X 10-^^ gcm-3 ( _1_ ] a(x) , (48) 

t 



where 



«(r, t) = 2.3 km/s ( j v{x) , (49) 

T \ { t 
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In Figs. 10 and 11, we plot the density and velocity profiles for the collapse from an initial 
state defined by A = 1.5 and T = 0.5, and a dynamics index 7 = 0.5. The solid curve in Fig. 

10 depicts the initial density profile, and the dotted curves represent the ensuing profiles 
(from top to bottom) at times t = 10^, 10^, and 10^ yrs. Likewise, the dotted curves in Fig. 

11 show the velocity profiles (from bottom to top) at times t = 10'^, lO"'', and 10® yrs. These 
results clearly illustrate the inside-out nature of the collapse, with gas inside the transition 
boundary 

/ t 

rBfti2.4pc — — , (51) 



IMyr J 

falling inwardly away from the overlying gas layers. 

From the onset of collapse, it takes a time ~ 10^ yrs for the formation of a young 
stellar object to occur, and an additional ~ 5 x 10^ yrs before stellar outflows disperse 
the surrounding core material. Thus, our model predicts that cores associated with young 
stellar objects would be characterized by densities ~ 5 x 10^ cm~^ and radii ~ 0.1 pc, in 
good agreement with observations (Jijina et al. 1999). In addition, the mass accretion rate 
per unit length is 

o \ 1/2 



M = IC\^—j mo = 1.2 X lO'^MoMyr-^pc-^mo. (52) 

An elongated core with a length of I — 0.5 pc would accrete a mass of Mace ~ lOM© in a 
time for a slightly overdense initial state (since mo ~ 0.15 in this case, as can be from 
the solid curve in Fig. 9), and Mace ~ 37Mq for A = 1.5 (mo = 0.6). These results are 
consistent with the low star formation efficiencies of 2 — 10 % deduced from observations, 
but seem to favor initial states that are only shghtly overdense. 



6. SMOOTH SINGULAR SOLUTIONS 

6.1. General Formalism 

The solutions presented in §3 did not cross through the singular surface, defined in the 
three- variable space x, a, and v as the surface on which T> = 0. While several different 
forms of solutions can pass through the singular surface (sec, e.g., Lou & Cao 2007), we 
focus here on those solutions which pass through it smoothly. That is, we consider solutions 
for which the derivatives v' — dv/dx and a' — da/dx exist on the singular surface. A 
necessary but not sufficient condition for the existence of these solutions is that V (and 
hence .4) also become zero, which occurs along a unique curve (referred to as the critical 
curve) on the singular surface for each set of values (a, 6*1,7). Alternatively, for solutions 
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which represent the collapse from an initially static, overdense state for which p — , a 
unique critical curve exists for each set of (A, F, 7) values. As a matter of illustration, we 
plot the (f c, and (cic, aj^) projections of the critical curve in Figure 12 for the parameters 
A = 1.5, F = 0.5, 7 = 0.5, with the upper panel presenting the reduced density profile, 
and the lower panel presenting the reduced velocity profile (dotted curves). The solid curves 
represent the solution for the corresponding collapse dynamics, which clearly does not cross 
the critical curve (or the singular surface). 

In order to properly treat solutions which pass through the singular surface smoothly, 
we obtain analytical solutions through a Taylor series expansion about the critical curve of 
the form x = Xc + 5, v = + and a = ac + where 5 « 1, and Vi and ai are 
the first derivatives in v and a evaluated on the critical curve. Since the system ODE's are 
second order, only the first two terms in the expansion are required. Values for vi and cti 
are obtained by substituting these terms into equations (17) and expanding in 5. Doing so 
then yields the relations 



and 



where 



-(!+-/>! + 



4F-l + 2(l-7)- 



^^1 + C„ - , (54) 



l-FV XcJ V 1/ xl 7 

The first derivative V\ is then the real, positive solution to the above quadratic equation. 
Full solutions can then be obtained by numerically integrating inward {5 < 0) and outward 
{5 > 0) from the corresponding Taylor series solutions. 

As an example, we consider the collapse from an initial state defined by the parameters 
A = 1.08, F = 0.25 and 7 = 0.5. Numerically integrating Eqs. (17) inward from an initial 
value of xq — 10^ (as described in §4), one finds that the equations become singular as 
X — > 0.336, i.e., as the solutions approach the singular surface. Since a real and positive 
solution vi exists for Eq. (54) at Xc = 0.336 for this case, this solution can thus pass 
smoothly through the singular surface (by crossing through the critical curve), and can 
be matched to the Taylor scries solution on the other side of the critical curve. Further 
numerical integration shows that this solution also crosses smoothly through the singular 
surface at a second point Xc — 0.0322, as illustrated in Figure 13. The sohd curve in the 
upper (lower) panel presents the reduced density (velocity) profile obtained for our solution, 
and the dotted curves represent the (etc, Xc) and {vc, Xc) projections of the critical curve for 
the case being considered. It is interesting to note that a approaches etc as x ^ 0. The 
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value of mo was also determined and is depicted in Fig. 9 by the solid square data point, 
and is consistent with the extrapolation of the dot-dash curve corresponding to the same (F, 
7) values. 



6.2. The 7 = 1 case 



As noted in §4, density profiles that are out of hydrostatic equilibrium by being every- 
where overdense (i.e., are of the form ai — Aa^, vi — 0) do not exist for the 7 = 1 case, 
regardless of the form for the initial equation of state as defined by a and Ci. This result is 
somewhat surprising given the apparent insensitivity that the solutions presented in §4 have 
on the choice of 7, but is borne out through numerical investigation. However, a different 
class of collapse solutions which passes smoothly through the singular surface docs exist for 
this case, and is presented here for completeness. A full analysis of all solutions with 7 = 1 
(e.g., wind solutions, shock solutions) is outside the scope of this work and will be presented 
elsewhere. 

To keep our analysis as general as possible, we describe the gas in terms of a and Ci 
throughout this subsection. It is easy to show that when 7 = 1, the critical curve on the 
singular surface becomes a straight line defined by the linear relation Vc — kxc and a constant 
reduced density 

(2 + 2a)(a + A;) 



where /c is a (real) solution to the quadratic equation 



- 1 



3 + 



4a 



k + 



4a + 2a^ + 



2a^ 



0, 



(56) 



(57) 



derived by imposing the condition that V — V — A — 0. Given our focus on collapse 
solutions, we plot the negative roots of Eq. (57) in Figure 14 as a function of a for values of 
Ci = 1.2,1.5,1.8,2.1 and 2.4. 

Solutions that pass through the critical curve can be obtained by numerically integrating 
inward and outward from both sides, using the analytical solutions obtained via the Taylor 
series expansion presented in §6.1 (where Vc — kxc)- We note that the values of vi and ai 
obtained by solving Eqs. (43) - (55) are now independent of - a consequence of the ODE's 
which describe the gas dynamics (Eqs. 17) being invariant to the scaling transformation 
X rjx, v^r]v,a^a,p^ rfp^ and m rfm for the case 7 = 1. We therefore set 
= 1 in our analysis without loss of generality. 



Through numerical exploration, we find a class of solutions that have the limiting form 
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OL oc x^/" and v oc x^^^"^^"" in the limit x oo, and a oc x"^/*^^ and f oc a;(2-Ci)/Ci the 
hmit X 0, where 1 < Ci < 2. To illustrate this type of collapse flow, we plot solutions 
for the case a = —1.05 and Ci = 1.5 in Figure 15. The upper panel presents the reduced 
density profile, and the lower panel presents the reduced velocity profile. The dotted curves 
in these panels represent the (etc, Xc) and (w c, Xc) projections of the critical line. Solutions 
for the real density and real velocity, scaled so that x = 1 corresponds to t = 1 Myr and r 
= 1 pc, are presented in Figs. 16 and 17 for times 0.01, 0.1, 1 and 10 Myrs. Clearly, this 
case represents a class of collapse solutions characterized by a nonzero initial inward velocity 
that slows with time at a fixed radius after the gas becomes isothermalQ. Indeed, the gas at 
large radii/early times {x —>■ oo) is moving inward with nearly constant velocity, though the 
magnitude of the velocity differs from one part of the flow to another (i.e., each "parcel" of 
gas moves with a nearly constant velocity, such that du/dt ^ 0, but the velocity field for the 
gas depends on r and t). The density increases as the gas falls inward (x — *• 0), leading to a 
rise in pressure that eventually becomes effective at slowing the inward collapse. The break 
in density at x ~ 1 in this case results from the slowing down of the infalling gas, rather 
than from an inward-out collapse associated with typical self-similar collapse solutions. This 
class of solutions could, in part, idealize the dynamics of a cylindrical structure collapsing 
after the loss of magnetic support via ambipolar diffusion, where it becomes isothermal as 
turbulence dissipates during the collapse. 



7. CONCLUSION 

The main goal of our work was to obtain self-similar solutions which describe the collapse 
of an initially stationary, cylindrically symmetric gas that is overdense from hydrostatic 
equilibrium by a factor A. Hydrostatic equihbrium profiles are easily found for an assumed 
static equation of state p = a'", as specified through the choice of the index F. Unlike 
previous works, we allow the equation of state to evolve during the collapse, under the 
condition that entropy is conserved along a streamline. Doing so allows the equation of state 
to change from its initial form (as defined by F), to a different polytropic form as defined by 
the dynamic index 7. Physical solutions describing this type of collapse require < F < 1 
and 7 7^ 1. We present solutions for which the system equations do not become singular, as 
well as solutions which pass smoothly through the singular surface. 

Our solutions clearly exhibit the broken power law profiles typical of self-similar collapse 



^ This class of solutions is analogous to the Type II solutions presented in Lou & Cao (2007) for a spherical 
geometry with 7 = 4/3 
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flows. We find that the spectral indices which characterize the reduced density solutions are 
not sensitive to the choice of 7, and hence have similar forms when the similarity variable 
X = Af^r << 1. In contrast, these solutions depend sensitively on the value of F, both in 
terms of shape when a; >> 1 and in overall normalization. Similar behavior is also clearly 
observed for the reduced velocity solutions. The insensitivity of the collapse dynamics to the 
value of 7 results from the inside-out nature of the collapse. Specifically, the break at a; ~ 1 
exhibited in the density and velocity profiles (as shown in Figs. 4-7) occurs as a result of the 
gas being in some part supported by pressure {x >> 1) evolving to a state approaching free- 
fall (x << 1) as a result of the loss of that pressure support. This break-point also denotes a 
transition from an initial (static) equation of state (p = a^) to a dynamic equation of state 
(p oc a"*), with the transition occurring within a narrow region around x ^ 1. In terms of 
real variables, this transition occurs at a boundary which moves outward through the gas 
as governed by the static equation of state. As such, the gas dynamics depend sensitively 
on the value of F (which, of course, also governs the density profile p(r) of the initial state). 
In contrast, the gas pressure within this boundary, as described by the dynamic equation of 
state (and in turn 7), has little effect on the overall collapse dynamics. 

Although self-similar collapse flows have no characteristic mass scale, the collapse flow 
feeds material onto the central star/disk system at a well-defined mass infall rate M. Since 
the infalling material always approaches free-fall conditions on the inside (in the limit r — > 
0), the reduced mass determines the size of the infall rate, and is therefore an important 
parameter in the collapse problem. Our analysis shows that, as expected, larger infall rates 
occur for initial states that are more overdense. In addition, stiffer static equations of state 
result in larger mass accretion rates. In contrast, the dynamic equation of state (as defined 
by 7) has little effect on the value of mo. 

Our results indicate that collapse from an overdense state initially at rest cannot occur 
if 7 = 1, regardless of the form for the initial equation of state. This result is somewhat 
surprising given the apparent insensitivity to 7 exhibited by all other solutions we obtain. 
However, we do find a different class of collapse solutions which passes smoothly through 
the singular surface when 7 = 1. These solutions are analogous to the Type II solutions 
presented in Lou & Cao (2007) for a spherical geometry with 7 = 4/3. 
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Fig. 1. — Self-similar solutions for the reduced density describing the collapse of a gas defined 
by the parameters a = —1.25 and 7 = 0.25. All solutions are integrated inward from an 
initial set of values Xj = 5 x 10'^ and a{xi) = 5.78 x 10"'^, and the five values of Vi as 
discussed in the text. All five solutions are within 0.1 percent of each other, illustrating that 
the self-similar density profiles are not sensitive to the initial values of the reduced velocity. 
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Fig. 2. — Self-similar solutions for the reduced velocity describing the collapse of a gas 
defined by the parameters a — —1.25 and 7 = 0.25. All solutions are integrated inward from 
an initial set of values Xi — 5 x 10^ and a{xi) — 5.78 x 10~^, and the five values of Vi as 
discussed in the text. 
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Fig. 3. — Log p - log a profiles for the solutions presented in Figs. (1) and (2), illustrating 
the transition during the collapse from a static equation of state (p = q;"+^ = q;°-^^) as shown 
in the lower left corner to a dynamic equation of state (p oc a''' oc q;°'^^) as shown in the 
upper right corner. 
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Fig. 4. — Self-similar solutions for the reduced density describing the collapse of an initially 
static gas overdense from its hydrostatic equihbrium state by a factor of A = 1.5. The four 
profiles correspond to the four different choices of static and dynamics equations of state, as 
defined by the specified values of 7 and F. 
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Fig. 5. — Self-similar solutions for the reduced velocity describing the collapse of an initially 
static gas overdense from its hydrostatic equihbrium state by a factor of A = 1.5. The four 
profiles correspond to the four different choices of static and dynamics equations of state, as 
defined by the specified values of 7 and F. 
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Fig. 6. — Same as Figure 4, but for a different set of values for 7 and F. We note that the 
soUd curve and short-dashed curve lie over each other in this Figure, and are therefore hard 
to distinguish. 
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Fig. 7. — Same as Figure 5, but for a different set of values for 7 and F. 
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Fig. 8. — The value of the reduced mass at the origin mo as a function of the static equation 
of state index F for A = 1.5 and three different values of the dynamic equation of state index 

7- 
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Fig. 9. — The value of the reduced mass at the origin mo as a function of the overdensity 
parameter A = 1.5 for the specified values of the equation and state indices F and 7. The sohd 
square data point corresponds to the results obtained for a solution which passes smoothly 
through the singular surface, as discussed in §5.1. 
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Fig. 10. — Density profiles illustrating the inside-out collapse of a cylindrical cloud initially 
at rest and overdense by a factor A = 1.5, and for which F = 7 = 0.5. The sohd curve 
represents the initial {t — 0) density profile p — ApE, and the dotted curves show the profiles 
(from top to bottom) at times t — 10^, t — 10^ and t — 10^ years. 
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Fig. 11. — Velocity profiles illustrating the inside-out collapse of a cylindrical cloud initially 
at rest and overdense by a factor A = 1.5, and for which F = 7 = 0.5. The dotted curves 
show the profiles (from bottom to top) at times t = 10^, t = 10^ and t = 10^ years. 
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Fig. 12. — The dotted curves denote the density (upper panel) and velocity (lower panel) 
projections of the critical curve for the parameters F = 7 = 0.5 and A = 1.5. The sohd curve 
in each panel denotes the corresponding self-similar solution, also shown in Figs. 4 and 5. 
This solution does not pass through the critical curve. 
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Fig. 13. — The dotted curves denote the density (upper panel) and velocity (lower panel) 
projections of the critical curve for the parameters F = 0.25, 7 = 0.5, and A = 1.08. The solid 
curve in each panel denotes the corresponding self-similar solution, which passes through the 
critical curve at Xc — 0.336 and at Xc — 0.0322. 
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Fig. 14. — The negative (real) roots of Eq. xxxxxx as a function of a for values of Ci = 
1.2, 1.5, 1.8, 2.1 and 2.4, as denoted in the figure. 
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Fig. 15. — The dotted curves denote the density (upper panel) and velocity (lower panel) 
projections of the critical hne for the parameters a — —1.05, Ci = 1.5. The solid curve in 
each panel denotes the corresponding self-similar solution which passes through the critical 
line at the point Xc — 1- 
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Fig. 16. — Density profiles associated with the self-similar solutions presented in Fig. 15, 
scaled by setting x^, = 1 when r = 1 pc and t = 1 Myr. The solid curves show the profiles 
(from top to bottom) at times t = 0.01, 0.1, 1 and 10 Myrs. 
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Fig. 17. — Velocity profiles associated with the self-similar solutions presented in Fig. 15, 
scaled by setting x^, = 1 when r = 1 pc and t = 1 Myr. The solid curves show the profiles 
(from top to bottom) at times t = 0.01, 0.1, 1 and 10 Myrs. 



